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Abstract 

To investigate the present situation of the merging in the southern outer region 
of Abell 85, we carried out long (~ 100 ks) observations with Suzaku, and produced 
an X-ray hardness ratio map. We found a high hardness ratio peak in the east side 
of a subcluster located in the south of the cluster; an X-ray spectrum of the region 
including this peak indicates a high temperature of ~ 8.5 keV. This hot spot has not 
been reported so far. We consider that this hot spot is a postshock region produced 
by the infall of the subcluster from the southwest. By using the Rankine-Hugoniot 
jump conditions for shocks, the Mach number and the infall velocity of the subcluster 
are obtained as 1.5 ±0.2 and 1950^280 km s -1 , respectively, in the case of merging with 
the subcluster from the southwest direction. By using the redshift difference between 
the A 85 and the subcluster obtained from optical observations, the angle between 
the line of sight and the direction of the motion of the subcluster is estimated to be 
75l£ de grees. We estimate the kinetic energy of the subcluster and the energy used 
for intracluster medium (ICM) heating to be ~ 10 63 and £8 x 10 60 erg, respectively. 
This shows that the deceleration of the subcluster by ICM heating has been negligibly 
small. 

Key words: galaxies: clusters: individual (Abell 85) — galaxies: intergalactic 
medium — shock waves — X-rays: galaxies: clusters 
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1. Introduction 



Clusters of galaxies are the largest objects in the universe, and consist of several tens 
or more galaxies that are bound together by gravity (Abell et al. 1989). They consist of 
galaxies, an intracluster medium (ICM) which is a high-temperature plasma gas, and dark 
matter. Numerical simulations show that many clusters of galaxies grow by merging during 
the formation of large-scale structures (e.g., Davis et al. 1985; Ricker 1998). Investigations 
of merging clusters are necessary not only to understand the evolution of clusters of galaxies, 
themselves, but also to tracing the structure formation of the universe. 

Merging clusters that violently come into collision are well studied because they exhibit 
interesting features, such as extreme variations in the temperature and the density distributions 
of ICM (e.g., Markevitch et al. 2002; Markevitch et al. 2004). However, such intensely merging 
clusters are very few, and mostly giant cluster cannibalize small clusters or groups of galaxies. 
Therefore, it is important to investigate such a general type merging clusters. Drastic heating 
of ICM by shock wave is not expected in such cases, and the early stage of merging occurs in the 
outer region with low X-ray luminosity. Therefore, highly accurate temperature determinations 
and photon-collecting power are required for investigating of such merging clusters. 

As a result of recent advances in X-ray observations, it is possible to find local high- 
temperature ICM regions that reveal the existence of merging in galaxy clusters (e.g., Govoni 
et al. 2004; Bourdin & Mazzotta 2008). However, even more accurate X-ray observations are 
required to investigate in detail the shock heating by merging. The X-ray imaging spectrometer 
(XIS: Koyama et al. 2007) and X-ray telescope on board the X-ray observatory Suzaku have a 
low background in the high-energy band and a large collecting area, respectively. Suzaku is cur- 
rently the best X-ray observatory for observing general merging clusters because its instruments 
are very effective in detecting hot, diffuse regions ,such as shocked regions. 

The Abell 85 cluster (hereafter, A 85) is a nearby (z = 0.055; Oegerle & Hill 2001) 
typical merging cluster in which a giant cluster is cannibalizing small clusters. A 85 has been 
observed so far by many X-ray observation satellites (e.g., ROSAT, Pislar et al. 1997; ASCA, 
Markevitch et al. 1998; Chandra, Kempner et al. 2002; XMM-Newton, Durret et al. 2003, 
2005, hereafter D05) and in other wavelength regions (e.g., Durret et al. 1998, Boue et al. 
2008). X-ray observations have revealed a bright subcluster located about 9' south of A 85. 
Figure 3 in D05 showed that a high-temperature region was found between the main part of 
A 85 and the subcluster, and they argued that this "impact region" might be thermalized by 
merging with the subcluster. However, this map did not cover the wide surrounding regions of 
the subcluster because of low signal-to-noise ratio (S/N). 

Based on the high statistics temperature map obtained from Suzaku observations, we 
investigated the merging situations and the motion of the subcluster. The Suzaku observations 
and data reduction are described in §2, and the results of data analysis are given in §3. We 
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discuss for the impact direction, velocity and motion of the subcluster in §4. Finally, the 
findings of the present study are summarized in §5. Throughout this paper, we assume that 
H = 70 km s" 1 Mpc" 1 , fi A = 0.73, fi M = 0.27, and q = 0.5. At the cluster redshift, z = 0.055, 
1' corresponds to 64.2 kpc. All quoted errors are at the 90% confidence level unless otherwise 
stated. 

2. Observations and data reduction 

2.1. Observations 

We observed A 85 located at (RA, Dec) = (00 h 41 m 37f8, -09°20'33") (Abell et al. 1989) in 
equatorial coordinates using Suzaku (Mitsuda et al. 2007) on January 2007 with a net exposure 
time of 99.2 ks. The Suzaku pointing position is ~ 8' south from the center (the peak of X- 
ray surface brightness distribution) of A 85, which is (RA, Dec) = (00 h 41 m 52?0, -09°25'43") 
in equatorial coordinates. This position has been termed the "impact region" in D05. The 
XIS field of view (FOV) is presented at the left side of figure 1 on the X-ray image taken by 
XMM-Newton. 

We obtained the XIS and the hard X-ray detector (Takahashi et al. 2007; Kokubun et al. 
2007) data, but the XIS data are used only in this paper because we perform a spatially resolved 
spectral analysis. The XIS is an X-ray CCD camera that consists of one back-illuminated (BI) 
sensor (XIS1) and three front-illuminated (FI) sensors (XIS0, XIS2, and XIS3). Each CCD 
covers a field of view of 18' x 18', and has a spatial resolution of about 2' in half power diameter 
(HPD). The BI sensor has a higher quantum efficiency in the soft energy band, whereas the FI 
sensors have lower instrumental non X-ray background (NXB) in the hard energy band. 

In this observation, the XIS2 data are not available because of the XIS2 anomaly 1 . The 
XIS was operated in the normal clocking and 3 x 3 or 5 x 5 editing modes, and the "spaced-row 
charge injection" (SCI: Uchiyama et al. 2007, 2009) option was applied. Energy resolutions 
(full width at half maximum; FWHM) of FI and BI data in this observation are 5.9 keV of ~ 
150 eV and ~ 170 eV, respectively. 

2.2. Data reduction 

We used reprocessed, unfiltered XIS event files (version 2.0) using HEASoft v6.5.1 
(Suzaku software vlO.0) and the calibration database (CALDB) released on 2008 September 
5. The energy correction was performed by using "xispi", and events were selected by using 
"GRADE=0:0 2:4 6:6" and "STATUS=0:65535". Good time intervals (GTI) were determined 
by the standard criteria referred to "Updated Gain Calibration for SCI-on Data Nov 15, 2007" 
2 . To remove hot /flickering pixels, we applied "cleansis" to the event file. The X-ray image 

1 http://www.astro.isas.ac.jp/suzaku/doc/suzakumcmo/suzakumcmo-2007-08.pdf 

2 http://heasarc.gsfc.nasa.gov/docs/suzaku/analysis/sci^ain_update.html 
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generated from the cleaned event file is shown in the right panel of figure 1. 

For spectral analysis, redistribution matrix files (RMFs) and ancillary response files 
(ARFs) were made by using "xisrmfgen" and "xissimarfgen" (Ishisaki et al. 2007). The XMM 
MOS1 image of A 85 was used as the source distribution, after subtracting the background, 
but not correcting the vignetting effect. Using a different source image, we confirm that the 
vignetting correction is insignificant in spectral analysis compared to the statistical errors on the 
results. The background components, the NXB and the extragalactic cosmic X-ray background 
(CXB), were estimated by the following procedure. The NXB images and spectra were made by 
using "xisnxbgen" (Tawa et al. 2007). For the CXB spectra, we referred to a method described 
in Sato et al. (2007). The ARF of a one-degree-radius uniform sky and the RMF were used 
to construct a CXB spectrum with the "fakeit" command of XSPEC (version 12.4.0ad). The 
spectrum model of the CXB was a power-law model with photon index T = 1.4 and Sx[2-10 
keV]= 5.97 x 10~ 8 erg cm -2 s" 1 sr" 1 , where Sx[2-10 keV] is the X-ray surface brightness in the 
2-10 keV energy band (Kushino et al. 2002), and iV H = 2.8 x 10 20 cm -2 was adopted as the 
hydrogen absorption column density (Hartmann & Burton 1997). Using the same requirement, 
the CXB events were created by "xissim" command to make the CXB images. 

A bright point source was detected at the position of (RA, Dec) = (00 h 41 m 57?5, 
-09°24'40") in the XMM image. We call this source DFL2005-17, since it is named "source 
17" in table 1 of D05. The contribution of this point source in spectral analysis is insignificant 
because the flux in the 2-10 keV energy band obtained from XMM-Newton data is less than 
10% of the flux in Suzaku HPD, indicated by the black circle in the right panel of figure 1. 
Therefore, we neglect its contribution for spectral analysis hereafter. Based on the ROSAT 
all-sky survey image (Snowden et al. 1995, 1997), the surface brightness of the Galactic soft 
X-ray diffuse background around A 85 is relatively low. Therefore, we apply a simple model 
for estimating this component. 

3. Analysis and Results 

3.1. Hardness ratio map 

In this subsection, we generate an X-ray hardness ratio map and examine the overall 
temperature structure of the cluster. The two FI CCD data (XIS0 and XIS3) are added together 
to increase the signal-to-noise ratio. The BI data was not used in this analysis because of their 
higher backgrounds at high energy above 5 keV. After subtracting the background components 
(NXB and the power-low CXB), we made hard (3-10 keV) and soft (0.5-2 keV) X-ray images. 
Figure 2 shows the obtained hardness ratio and its relative error maps. 

Between the main cluster and the subcluster, there is a hard region that was reported 
previously (D05). Moreover, we found a much harder spot on the east side of the subcluster. 
Note that this is ~ 2' away from the point source, DFL2005-17. Hereafter, we refer to this spot 
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Counts Counts 

Fig. 1. X-ray images of A 85. Left: Image obtained with XMM-Newton/MOSl in the 0.5-10 kcV 
band with background subtraction, but without any vignetting correction. Smoothing over 30" has been 
applied. Contour levels are logarithmically spaced. The magenta square represents the field of view of 
the Suzaku/XIS. The orange dashed ellipse indicates the "impact region" defined in D05. The white cross 
indicates the position of the X-ray point source "DLF2005-17" . Rright: The Suzaku/XISO image in the 
0.5-10 keV energy band. Smoothing over 17" has been applied. The white cross is the same as that for the 
left panel. The black circle indicates the size of the Suzaku point spread function (half power diameter of 
2'). No background has been subtracted and no vignetting correction has been performed. The calibration 
source regions are excluded. 

as the hot spot. Although a part of this hot spot may have been seen in the XMM-Newton 
results, main parts of the spot are not included in the analysis region of D05 (see figure 3 
in D05). The relative error map indicates that the detection of the hot spot has statistical 
significance (with a relative error < 15%). Therefore, the Suzaku observation detected this hot 
spot robustly for the first time. 

3.2. Spectral analysis 

3.2.1. Temperature of each region 

We divided the field of view into seven regions based on the hardness ratio map, as 
shown in the right panel of figure 3. We extracted the X-ray spectrum from these seven regions 
separately. 

The energy ranges used for the spectral analysis were 0.5-10 keV for FI and 0.5-7 keV 
for BI. The energy range of the anomalous response around the Si K-edge (1.825-1.840 keV) has 
been excluded in both cases. To describe the cluster emission, we used a collisional ionization 
equilibrium model (APEC; Smith et al. 2001) modified by the Galactic photoelectric absorption 
(wabs in the XSPEC), using XSPEC. The Galactic hydrogen column density was fixed to the 
reported value, 2.80 xlO 20 cm~ 2 (Hartmann & Burton 1997). To model the Galactic soft X- 
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Fig. 2. Hardness ratio (3-10 keV/0.5-2 keV; left) and relative error (right) maps of the FI. The pixel 
width is 1' and smoothing over 1' has been applied. The contours and the white cross are the same as 
those in figure 1. The background is subtracted without any vignetting correction. 

ray background, another APEC model was added to the model. By fitting the spectrum in 
"Blank" region, temperature, kT, is determined as 0.18 keV. This value was used for other 
regions. In contrast, normalizations of other regions are left free to compensate the variation 
of the background flux. 

Figure 4 shows the fitting result of each region and table 1 summarizes the obtained 
parameters. The fitting from the BI gave larger statistical errors than those from FI, in general. 
The spectral fitting of "Main" region indicates an excess at an energy of around 1.1 keV, as 
shown in upper-left panel in figure 4. Adding another thermal component, we could improved 
the fit. This is possibly because of the cool component commonly seen in the center of galaxy 
clusters (e.g. Fabian et al. 1984). 

The left panel of figure 3 shows the temperature profiles obtained by using the three 
detectors. The temperature of "East" region from the XIS3 is inconsistently higher than that 
obtained from other detectors. This is possibly caused by a calibration error at the edge of 
CCDs. We omit any discussion of "East" region hereafter. The "Middle" and "Hot" regions are 
significantly hotter (kT ^7 keV) than the "Main" and "Subcluster" regions (kT ~6 keV). This 
is consistent with the hardness ratio map (figure 2). A more detailed temperature distribution 
is given in the next subsection. 
3.2.2. Detailed temperature map 

In order to investigate the detailed temperature distribution of the cluster, we made a 
fine temperature map focusing on the hot component in particular. After excluding the "Blank" 
region with low S/N, we divided the whole region into 2' x 2' boxes. We perform spectral fitting 
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Fig. 3. Left: Temperature profiles of the XISO (black), XIS1 (red), and X1S3 (green) detectors. The 
horizontal axis represents regions extracted spectra. Right: The spectral extraction regions are indicated 
on the XISO image. The contours of the hardness ratio map (the left panel of figure 2) are also overlaid. 
The "Radio" region includes some radio sources (Bagchi et al. 1998; Lima Neto et al. 2001). The white 
cross is the same as that in figure 1. 



Table 1. Results of spectral fitting for the seven regions, using the added FI (XISO + XIS3) spectn 



Region 


Main 


Middle 


Hot 


Subcluster 


Radio 


East 


Blank 


fcT(keV) 


c 70+O.O6 


7 ao+0.32 
' ' zo -0.31 


6.87l°;« 


5 76+ - 22 

°-' D -0.21 


r ro+0.30 

°- oo -0.28 


6 7^+ - 50 

°-'°-0.32 


5 07+ - 18 


Redshift 


fl 562+ - 0005 

U.UOD/_ 0005 


U.UD4U_ 0025 


0.056218:582 


0581 + ' 0036 
u.uooi_ 0026 


0565+ - 0048 
u.uooo_ 0046 


U.UODJ_ 0027 


6 03+ - 0060 
u.uouo_ 0051 


zT(solar) 


46+ - 02 


n o7 +0.05 
u - z '-0.05 


26+ - 06 

u - zo -0.06 


21 + - 04 

U - Zi -0.04 


26+ - 06 

u - ZD -0.06 


16+ 010 


u - i '-0.04 


Norm*'''' 


38 - -cu 


9 9 r;+0.03 
z - ZJ -0.03 


1 -?n +0 - 03 

i - ou -0.03 


9 44+0.04 

z - q4 -0.04 


i C9+0.04 
i - dz -0.03 


i 17+0.02 
1,A 1 -0.03 


9 OO+0.05 

z -°°-0.05 


rVd.o.f. 


1368/938 


356/303 


416/362 


309/276 


339/311 


234/249 


305/303 



* The APEC model normalization is 10" 17 § n e n H dV / (birD A {1 + z) 2 )- 

t Normalizations are rescaled multiplying by SOURCE_RATIO_REG in FITS header of the calculated ARFs using "xissi- 
marfgen" . 



in the similar way as discribed the previous subsection, but only used the FI spectra. 

The half power diameter of the telescope point-spread function, 2', is comparable to the 
spatial bin size in our analysis. We estimated that about 50% of the event in a box originated 
from surrounding regions (see also Ota et al. 2007). Therefore, our obtained temperature map 
was substantively diluted, and the actual temperature variation could be larger than the present 
result. The reduced x 2 of spectral fittings in all boxes are in the range of 0.86-1.20, except for 
CCD edge regions (boxes 00-50). 

The obtained temperatures are shown in figure 5. The cluster center (box 03) and the 
subcluster (boxes 44, 45, 54 and 55) are cool {kT ~ 4.8 keV and ~ 4.9-6.8 keV, respectively) 
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relative to a high-temperature spot between the main part of the cluster and subcluster, (box 
34, kT ~ 8.2 keV, hereafter we refer to this spot as the middle spot) and the hot spot (boxes 
43, kT ~ 8.5 keV). The low temperatures at the cluster center and the outer envelope ("Blank" 
region) may reflect the general temperature profile of galaxy clusters (e.g. Allen & Fabian 
1994). This temperature distribution is consistent with the hardness ratio map obtained in 
§3.1. 

The same results are shown in figure 6 in one-dimensional plots along the detectors' 
axes. The temperature profile of "boxes 40-45" in the figure shows a significant temperature 
difference of ~ 1.7 keV between boxes 43 (hot spot) and 44 (a part of the subcluster). Other 
temperature difference can be seen between boxes 03 and 04, between boxes 13 and 14, and 
between boxes 53 and 54. Some of the temperature variations are already shown in figure 3 
of D05, except for the difference between boxes 53 and 54. This analysis have revealed the 
temperature structure including the hot spot very robustly, and in wider spatial area than in 
previous measurements. 
3.2.3. Maps of other quantities 

In search for further evidence of the merging of subsystems and features in A 85, we 
produced projected maps of particle number density, pressure, and entropy index, as shown 
in figure 7. The electron number density is calculated from the normalization of the APEC, 
K = 10~ 14 / nnn e dV/(4:7iD\(l + z) 2 )(cm" 5 ), where n e and rin are, respectively, the electron and 
hydrogen densities in units of cm -3 and Da is the angular diameter distance in units of cm. 
The solid angle AQ of each box is 4 arcmin 2 . Then, projected n e is obtained as (Henry et al. 
2004) 

n e = 7.28 x 10~ 7 (1 + z) ( -^-) ( ?—^) ^ ( — ^ ^ cm" 3 . (1) 

v ; \lMpcy V10- 3 cnW yiMpc 3 / v; 

Here, we assumed that the volume is V ~ (2/3) D\ AQ(R^ ir — R 2 ) 1 ^ 2 , where R is the projection 
radius from the center of the cluster and R vn is the virial radius, ~2 Mpc (Reiprich & Bohringer 
2002). By using H:He = 9:1 in terms of numbers of atoms as the abundance ratio of the ICM, 
the total particle (electron plus ion) number density is obtained to be n = (13/ll//)n e (White & 
Silk 1980), where \i ~ 0.62 is the mean molecular weight. The pressure is given by P = nkT, and 
we use the entropy index S = kTn~ 2 ^ 3 instead of the thermodynamic definition (3fc/2)ln[Tp~ 2 / 3 ], 
where p g is the gas density. 

The number density map naturally resembles the X-ray surface brightness distribution. 
The high-temperature region (boxes 33, 34, 43 and 53) has a relatively high pressure. The east 
side of the subcluster (boxes 42, 43, 52 and 53) has a higher entropy than the surrounding 
region. The entropy in this region might be increased by merging. These maps of particle 
number density, pressure and entropy are used in the discussion below. 
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3.3. Complementary analysis using optical data 

Although we expected that a significant in difference redshift between the main part of 
A 85 and the subcluster would be detected in the X-ray observation, it was not obtained (see 
table 1). Then, we consider that the difference may be detected by using the radial velocities 
of the member galaxies obtained from the optical spectroscopic observations whose accuracy of 
the redshifts are higher than those derived from X-ray data. In addition, we may confirm other 
features of the member galaxies using such information as optical colors and spatial distribution. 

The catalog of Durret et al. (1998) is used for our analysis. This catalog includes the 
information on the positions and the radial velocities of member galaxies obtained from their 
optical spectroscopic observation. The catalog covers about 1 square degree of cluster center, 
and the typical error on radial velocities in this catalog is smaller than 100 km s" 1 . Though the 
catalog does not cover homogeneously the whole region, it is enough for identifying physical 
concentration of galaxies around main and subcluster regions, which is essential for deriving 
mean radial velocities of both components. Since the catalog includes not only early-type, but 
also late-type galaxies, we can perform the color analysis for the member galaxies by adding the 
color information to the galaxies. For the purpose, we matched the catalog with the photometric 
catalog of the Sloan Digital Sky Survey (SDSS; York et al. 2000) DR-7 (Abazajian et al. 2009) 
by identifying the SDSS galaxies within the position of ±2 arcsec. In this subsection, we use 
that data set hereafter. 

The spatial distribution of the member galaxies is shown in the left panel of figure 8. We 
considered that the member galaxies of A 85 should have radial velocities cz within 13500-19500 
km s _1 because the systemic velocity of A 85 is 16507 ± 102 km s" 1 (Oegerle & Hill 2001). To 
confirm the radial velocities of individual regions visually, we separate the member galaxies into 
three radial-velocity groups. Blue, green and red points indicates the range of the radial velocity 
of 13500-15500, 15500-17500 and 17500-19500 km s" 1 , respectively. Because of its compact 
concentration and positional coincidence, we here considered that the 10 galaxies shown in the 
right panel of figure 8 surrounded by dotted line are the member galaxies of subcluster. In the 
region between the center and the subcluster, there are some galaxies with relatively large radial 
velocities, shown by red points. It might be indicative that there was another collision by small 
group of galaxies previously in this region. Southeast of the subcluster, there is a clustering of 
galaxies with small radial velocities, which is consistent with the report that an X-ray filament 
expanding in this direction (Durret et al. 2003). The member galaxies of the subcluster (see 
the right panel of figure 8) listed in table 2 have almost the same radial velocities ,and their 
mean radial velocity is 17028 ±216 km s _1 . This value is lager than the systemic velocity of 
A 85 ~ 520 ± 240 km s _1 , and it is consistent with the fact that a significant difference radial 
velocity could not be detected by this X-ray observation. 

The diagram, u — g versus g, of the member galaxies is shown in figure 9. As discussed 
in Bower, Lucey, and Ellis (1992), there must be "red-sequence" made by early- type member 
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Table 2. Coordinates, magnitudes of the u and g bands, u — g, and cz with errors for the member galaxies of the subcluster. 



RA 


Dec 


u 


9 


u~g 


cz 


(J2000.0) 


(J2000.0) 








(km s _1 ) 


00 41 39.0 


-09 27 48 


19.99 


18.30 


1.69 


17257±52 


00 41 39.6 


-09 27 31 


18.58 


17.16 


1.42 


16647±62 


00 41 43.1 


-09 26 22 


17.37 


15.41 


1.96 


16886±35 


00 41 48.1 


-09 27 30 


19.49 


17.79 


1.70 


17203±66 


00 41 49.1 


-09 29 03 


19.78 


17.80 


1.98 


18437±136 


00 41 50.4 


-09 28 18 


19.93 


18.09 


1.84 


17293±49 


00 41 51.9 


-09 31 58 


19.45 


17.58 


1.87 


16523±243 


00 41 52.3 


-09 30 16 


17.51 


15.62 


1.89 


17164±33 


00 41 53.2 


-09 31 17 


19.56 


17.81 


1.75 


17121±67 


00 41 53.3 


-09 29 29 


18.01 


17.14 


0.87 


15751±75 



galaxies in rich clusters of galaxies. From the distribution of the member galaxies, we can derive 
the relation, u — g = — 0.061c/ + 2.9 (dotted lines), and it is very consistent with that of Coma 
cluster (see left bottom of figure 9). The facts shown in this figure are as follows: (1) From the 
photometrical data of SDSS-DR7 only, it is difficult to find the "red- sequence" . As there are 
some reports that some components are overlapping along the line-of-sight in the field of A 85 
(e.g. Bravo- Alfaro et al. 2009), the "red-sequence" might not clearly be seen. (2) Many of the 
member galaxies gather near the relation, while there exist a certain number of galaxies below 
them (bluer than them). This implies that A 85 is enriched with blue galaxies, which may be 
proceeding more actively with star formation than typical rich clusters. On the other hand, we 
could not get any remarkable features in the spatial distribution of the member galaxies with 
color information. 

4. Discussion 

4-1. Impact direction 

Previous studies suggest that the subcluster has moved from the southeast of the present 
position, based on the following findings: (1) X-ray observations show that the cold front is 
located between the subcluster and the main part of A 85 (Chandra, Kempner et al. 2002). 
(2) An X-ray filament extends from the subcluster to the southeast direction (XMM-Newton, 
Durret et al. 2003). (3) The "impact region" has a high temperature (XMM-Newton, Durret 
et al. 2005). (4) Some Ha emitting galaxies are in the southeast of the subcluster (Boue et al. 
2008). 

However, we found a hot spot in the east side of the subcluster. This feature can not 
be easily explained by merging of the subcluster from the southeast direction. Instead, it is 
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more plausible to consider that this hot spot is caused by merging from the southwest direction. 
Therefore, we proposed and examined a new model in which the subcluster collides from the 
southwest (see the left panel of figure 10). Regarded from this new point of view, the "impact 
region" appearing in D05 is considered to be the northwest end of the shocked region. Moreover, 
the listed facts wiht which the collision from the southeast direction is observationary supported 
do not critically conflict with merging from the southwest direction. 

4-2. Motion of the subcluster 

Here, we estimate the infall (merging) velocity of the subcluster using Rankine-Hugoniot 
jump conditions (RH conditions) in the case of merging with the subcluster from the southwest 
direction. To obtain the information about the motion of the subcluster, we consider the 
boxes 41 and 43 to be the typical parts of preshock (index 1) and postshock regions (index 2), 
respectively (see the right panel of figure 10). 

The way we consider this problem is as follows: (I) "The shock-wave speed is equal to 
the infall velocity of the subcluster." Since the shock wave is formed in the open system, we 
consider that the whole structure of the shock wave plus the subcluster keep a steady form, 
i.e. the shock-wave speed is equal to the present velocity of the subcluster relative to A 85. 

(II) "The postshock region is consists of two temperature components." The results obtained 
in previous sections are physical quantities integrated along the line of sight. In order to 
obtain deprojected quantities, we fit the two-temperature model ((APEC + APEC) x wabs + 
APECg*) to the spectrum in the postshock region (box 43). We consider that the quantities 
of high-temperature {kT 2 = 10.4jlj;| keV) component are those of the postshock region. The 
extent and the depth of the postshock region are estimated to be 4 arcmin 2 and 5-8 arcmin 
when the shock wave passed through the region. This gives the deprojected particle number 
density as n 2 = 1.70 ±0.20 x 10~ 3 . The estimations of the volume ot the postschok region 
are consistent with the emission integral of the region. (Ill) "The initial conditions of the 
postshock region are equal to those of the preshock region." Regions in the same radius would 
be similar physical state in relaxed galaxy clusters before the merging. Since the postshock 
and the preshock regions are seem to be at almost equal distance (~ 8') from the center of A 
85, we consider that initial conditions in these regions were nearly equal. 

Using the obtained temperature and the particle number density of the postshock region, 
we estimate the Mach number related with the merging, and the infall velocity of the subcluster. 
We later note that the infall velocity is somewhat affected by a possible breaking of assumption 

(III) as the difference in the initial conditions between the preshock and the postshock region 
before the merging. 

By using RH conditions, the following relation is obtained (Landau & Lifshitz 1959): 



P2 2 7 



Ml 



7-1 



(2) 



Pi 7 + 1 



7 + 1 
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where 7 = 5/3 is the adiabatic index for the fully ionized plasma and Aii — v\/c s \ is the 
Mach number, and where v 1 and c s i are the infall velocity of the subcluster and the sound 
speed in the preshock gas, respectively. P\ and P 2 are given by Pi = 1.09lo!i2 x 10~ n erg 
cm" 3 (kTx = 6.57±8;|? keV, n x = 1.04 ±0.03 x 10~ 3 cm" 3 ), and P 2 = 2.82±g;|f x 10" 11 erg cm" 3 
(kT 2 = lOAtH keV, n 2 = 1.70 ±0.20 x 10" 3 cm" 3 ). We then obtain P 2 /P 1 = 2.59^;|, c sl = 
1300^70 km s _1 ; Mi and v\ are estimated as 1.5 ±0.2 and 1950l2go km s^ 1 . By using the 
difference redshift between the subcluster and A 85, 520 ± 240 km s" 1 obtained in §3.3, the 
angle between the line of sight and the direction of motion of the subcluster is found to be 
cos _1 (520 ± 240/1950^280) — 75lg degrees. Here, we estimate the effects of a possible breaking 
of the assumption (III). From v\ = Ai±c s i, we found that V\ mainly depends on (P 2 /ni) l l 2 . 
Therefore, for example, V\ decreases (increases) by 5% if n\ increases (degreases) by 10%. 

Finally, we roughly estimate the present kinetic energy of the subcluster and the kinetic 
energy loss by the ICM heating and compare both. The former, K s , can be estimated to 
be (l/2)M s vf. Here, M s is the mass of the subcluster, which is obtained from the balance 
its self-gravity and pressure gradient, i.e. ~ kT s R s / fiGm-n ~ 3 x 1O 13 M , where T s is the ICM 
temperature, the R s radius of the subcluster, G the gravitational constant, and mn the hydrogen 
mass. Referring to the above discussions, K s is estimated to be ~ 10 63 erg. On the other 
hand, the energy loss by the ICM heating is estimated to be /S.E\ iCat = kAT^ e&t x V^ept X n, 
where &;ATi lcat is the difference temperature between the preshock and the postshock region, 
Kwept = ^-^s-^tr the total volume of the region swept by the subcluster, and n the mean particle 
number density of the region. D tT is the distance along which the subcluster has traveled. 
Using the values, A;AT hcat < 7.4 keV, R s = 150 kcp, D tT = 500 kpc and n = 10~ 3 cm" 3 , we obtain 
APhcat ^ 8 x 10 60 erg. The fact that APheat is very much less than K s means the subcluster is 
scarcely decelerated by ICM heating. 

5. Summary 

The high-statistics map of the X-ray hardness ratio obtained from long-time (~ 100 ks) 
observations of A 85 with Suzaku revealed a large peak on the east side of the subcluster. From 
the spectral fittings, a comparable result was obtained. A detailed temperature map shows a 
higher temperature (~ 8.5 keV) region that overlaps with this peak, and a similar feature can 
be seen in the maps of pressure and entropy. This high-temperature region (hot spot) has not 
been reported so far. From the location of the hot spot, we consider that a postshock region 
produced by the infall of the subcluster from the southwest, in contrast to previous studies, 
which have recognized that the subcluster is inf ailing from the southeast. 

By using the Rankine-Hugoniot jump conditions for shocks, Mach number and the infall 
velocity of the subcluster are estimated to be 1.5 ±0.2 and 1950^280 km s _1 , respectively. On 
the other hand, by using the redshifts of the subcluster's member galaxies obtained from optical 
observations, the radial (line-of-sight) velocity difference between the subcluster and A 85 was 
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found to be 520 ±240 km s , in excess for the subcluster. These velocities yield the conclusion 
that the angle between the line of sight and the direction of the motion of the subcluster is 
75l£ de grees. The present kinetic energy of the subcluster and the energy loss by the ICM 
heating are roughly estimated to be ~ 10 63 and ^ 8 x 10 60 erg, respectively. This means that 
the subcluster is scarcely decelerated by ICM heating. 

We thank an anonymous referee for useful comments. Data analysis was carried out in 
part on the common-use data analysis computer system operated the Astronomy Data Center 
(ADC) of the National Astronomical Observatory of Japan. S.M. acknowledges support from a 
Kyoto Sangyo University Research Grant. We thank all the Suzaku team members for operating 
the satellite and assistance with data analysis. 
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Fig. 4. Results of spectral fitting. Each panel shows the background (the power-low CXB + NXB) 
subtracted spectra of the region indicated in the right panel of figure 3. Black, red, and green crosses 
indicate XISO, XIS1, and XIS3 data, respectively. In the upper part of panel, the solid lines indicate the 
best-fit model for "APEC x wabs + APECg" , and the lower dotted lines show the component "APECg" 
of the model, as the Galactic soft X-ray component. Note that in "Main"- and "East"- region fitting, the 
components do not appeared in each plot because of the low flux. In the lower part of each panel, the fit 





kT[keV] 

Fig. 5. Temperature distribution map. The contours of X-ray surface brightness from the Suzaku data 
are shown by black lines on a logarithmic scale. The contours of hardness ratio are indicated by red dashed 
lines on a liner scale. Calibration sources located at northwest and southwest have been excluded. The 
white cross is the same as that in figure 1. The box address map is given in the bottom- left, where "CS" 
represents the positions of the calibration source regions. 
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Fig. 6. One-dimensional temperature profiles along the east-west direction. The box addresses are shown 
in figure 5. 
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Fig. 7. Particle number density, pressure and entropy maps are shown in top left, top right and bottom 
left, respectively. The contours are the same as those in figure 5, and the white cross is the same as that 
in figure 1. 
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Fig. 8. Left: Spatial distribution of member galaxies with radial velocity. The contour levels are the 
same as those in the left panel of figure 1. The part enclosed by dashed lines is enlarged in the right panel. 
Right: Close-up of the left panel. The points and the contours are the same as those in the left panel. It 
is conjectured that galaxies enclosed by red dashed lines are the member galaxies of the subcluster, listed 
in table 2. 
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Fig. 9. Color-magnitude diagram, u~ g versus g, of A 85 field. The dotted line represents the relation 
u — g = — 0.06f g + 2.9. Green crosses indicate galaxies distributed along the relation within the u — g color 
difference of ± 0.2 mag. Red and blue crosses represent galaxies that are redder and bluer than those 
shown by green points, respectively. Black points indicate the photometric data of SDSS DR7 within the 
field covered by Durret et al. (1998). A similar diagram for Coma cluster is shown in the small panel of 
the left bottom. In the small panel, points indicate photometric data of SDSS DR7 in Coma cluster, and 
the dotted line represents the relation u — g = —0.0649 + 2.8. 



20 



Counts 

4.0 5 S 7.1 M.6 10 

Fig. 10. Left: Impact directions of the subcluster proposed in previous studies (white arrow) and in this 
study (yellow arrow) superimposed on the XMM-Newton image. Right: Bird's eye view of the shock wave 
and the motion of the subcluster. The yellow arc indicates the expected location of the shock front and 
the ellipse indicates the subcluster. The temperature map is the same as that in figure 5. 
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